Preserving Entanglement of Flying Qubits in Optical Fibers by 

Dynamical Decoupling 
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We theoretically investigate the influence of dynamical decoupling sequence in preserving entan- 
glement of polarized photons in polarization-maintaining birefringent fibers(PMF) under a classic 
Gauss 1/f noise. We study the dynamic evolution of entanglement along the control sequence em- 
bedded fibers. Decoherence due to dispersion of polarization mode in PMF can be dramatically 
depressed, even for a wild optical width. Entanglement degree can be effectively preserved while 
the control pulse is implemented. 



Entanglement plays a central role in quantum commu- 
nication process. Due to an interaction with the uncon- 
trollable degree of freedom of the environment, quantum 
system may lose its entanglement degree. In an optical 
fiber-based quantum communication channel, the resid- 
ual optical birefringence randomly accumulating along 
the fiber set up a main obstacle to maintain the fidelity 
of information. Polarization-entangled photons pairs dis- 
tributed over optical fibers even suffer the process of 
abrupt disappearance of entanglement known as a 
phenomenon of entanglement sudden deathj^, Thus 
it is upper most important to find effective methods for 
preserving entanglement of polarized photons propagat- 
ing in optical fibers. 

Several approaches have been developed to address this 
issue. Notable examples are decoherence-free sbuspaceJU, 
5] and quantum error-correction codes [1,0]: both based 
on carefully encoding the quantum information into a 
wider, while partially redundant, Hilbert space. Alter- 
native approach is quantum feedback Q. In this tech- 
nology information channel is designed to be closed-loop 
with appropriate measurements and real-time correction 
to the system. However, all these strategies have the 
drawback of requiring a large amount of extra resources. 
On the contrary, an open-lope control method known as 
dynamical decoupling avoids all these hindrances. 

Dynamical Decoupling(DD) is a simple and effective 
method for coherence control. In this technology un- 
desired effects of the environment are eliminated via 
strong and rapid time-dependent pulses faster than the 
environment correction time. The physical idea behind 
DD scheme comes from refocusing techniques in Nu- 
clear Magnetic Resonance (NMR) systems [9( and then 
be extended to any other physical contexts in the last 
decades, such as nuclear-quadrupole qubits and electron 
spins qubits [ll| . Preserving of entanglement between 
two stationary qubits by dyn amical decoupling also has 
been wildly discussed [llHl3j- Some prominent exam- 
ples of DD schemes are the wildly used Hahn's spin 
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echo, periodic DD (PDD), Carr-Purcell-Meiboom-Gill 
(CPMG)Q, concatenated DD (CDD)[lJ], Uhrig DD 
(UDD) 16], and recently proposed near-optimal decou- 
pling (QDD)[l7| used to eliminate general decoherence 
of qubits. 

Extension of the time-dependent DD sequences to the 
spare-dependent dynamic evaluation process lead to the 
idea of applying DD controls into the optical fiber-based 
quantum communication channel [Lsl Il9j . Experimen- 
tal implementation has notably revealed the potential 
of this extension [2(| [2l|. However, this delicately de- 
signed experiment is limited to the condition that the 
noises are systematically introduced in a non-stochastic 
way, while in a real fiber-based channel they are dis- 
tributed randomly. In this paper, we investigate the 
performance of CPMG sequence in preserving entangle- 
ment of polarized photons under stochastic classic Gauss 
noises. We consider here the entangled photons are dis- 
tributed through polarization-maintaining(PM) birefrin- 
gent fibers, in which case the polarized photons suffer 
a pure dephasing process, Tl >> T2. Such building 
blocks of coherence control for a dephasing dynamic can 
be extended to general decoherence processes in ordi- 
nary single mode optical fibers with more complicated 
DD schemes like CDDpJ] or QDDpjj]. In the following, 
we show the derivation of entanglement evolution using 
the filter-design method [22| and give the numerical sim- 
ulation with some special initial states. 

We consider the situation of photon distribution de- 
picted in Fig. 1. Photon A is preserved by a quantum 
register, while the entangled photon B, carrying the en- 
coded information propagates through the fiber to the 
recipient. The fiber is embedded with 7r pules realized 
by half- wave plates. In our scheme the intervals between 
waveplates are fixed to a certain scale in the fiber and the 
waveplates sequence are arranged as continues CPMG cy- 
cles. Note that the CPMG sequence has a self-construct 
feature as illustrated in Fig.l. In other words we can 
view two cycles of CPMG sequences with N pules each 
as one cycle of CPMG sequence with 2N pules. This fea- 
ture of CPMG sequence allows us to analysis the dynamic 
evaluation of entanglement through the fiber directly, as 
illustrated in the following. 
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\w, — lj > is the frequency basis vector of photon pairs, 
w denotes the offset from the central frequency. \P > 
represents polarization modes, which can be expressed as 
the superimposing of horizontal and vertical polarization 
basises \HH >, \HV >, \VH >, \VV > 

Due to the stochastic dephasing noise described in 
Hamiltonian (1), the output state of the photon pairs af- 
ter photon A propagating through the fiber for a length 
L without DD sequence implemented is: 



2L 



FIG. 1: (a)Photon distribution scheme. (b)The length be- 
tween pulses embedded in the fiber is fixed to a certain scale. 
Measuremnet at length L in the fiber corresponds to a 2 pulses 
CPMG sequence, while measurement at length corresponds to 
a 4 pulses CPMG sequence. 



Under a photon distribution process described above, 
the qubits-environment interaction Hamiltonian can be 
written as L 23j: 

H = i J duj(\uj A >< w\ ® b(oj, L)Oz + \lo a >< u\ ® I B ) 

(1) 

In nondispersive media only the first order of u> remains 
in L)[23|. Thus we can rewrite b as: 



b(u,L) =u)(Sl + P(L)) 



(2) 



fHs a constant and f3(L) presents for the stochastic fluc- 
tuation of the noises with a zero mean, and it has a two- 
point correlation function: 



S(L 1 -L 2 )=<0(L 1 ),0{L 2 ) > 



(3) 



Where <> corresponds to the average with respect to 
the noise realizations. The statistical properties of the 
environment can also be expressed as the spectral density 
of noise: 



S(u) 



t S{L)dL 



(4) 



In the following discussing the statistics of fluctuations is 
assumed to be Gaussian, in which case the noise is com- 
pletely defined by the first-order correlation function S(L) 
in equation(3). We assume the polarization entangled 
photon pairs are generated via spontaneous parametric 
down conversion. The quantum state of the generated 
photon pairs can then be written as[l|: 
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(6) 



where the subscript s stands for a given realization of 
noise before taking the ensemble average. The den- 
sity matrix that characterizes the detected polarization 
qubits can be abtended by a process of tracing over the 
frequency modes and taking ensemble average of the en- 
vironment noise: 



p = J Df3{T 



r"\ipout >< VwD) 



(7) 



The integration is a Gaussian functional integral with the 
variable ft(L) having a Gaussian distribution. 

The elements of the resulting density matrix are then 
given by: 



Pii,12,34 = Pii, 12, 34(0), 

i = 1,2,3,4 



y/l + o*f(L) ^ l + &*f(L)' 
i,j = others, (8) 

where f(L) can be expressed as the integral of the noise 
spectrum S(w): 



(9) 



F(x)is a filter function depends on the shape of DD se- 
quence. For a free evaluation without control imple- 
mented F(x) = 2 sin 2 — . The fiber-based DD control can 
2 

be considered as a filter design process(See Ref.[22| for 
details). Any kinds of control sequence can be expressed 
as a particular filter function. Analytical expression for 
CPMG sequence with pulses of even number N(replace 
sin 2 x with cos 2 x for odd-N CPMG) embedded before 
the measurement point can be written as: 



Fix) =8sin 4 ( — ) s in 2 (-)/cos 2 — 



(10) 



Note that the intervals between pulses are fixed in our 
scheme. Namely, N is proportional to the length of the 



3 



0.35 
0.3 

CD 25 
O 

CD 0.2 



2 pulses in each unit length 
4 pluses in each unit length 
6 pluses in each unit length 
8 pluses in each unit length 



(a) 



40 50 60 

Length of fiber 




Length of fiber 



FIG. 2: (a)Evoluation of entanglement along the fiber for 
different pulse intervals. (b)Evoluation of entanglement with 
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FIG. 3: Concurrence as a function of the total pulse number 
within 50 units length. 



fiber. As the particular property of CPMG described 
above, we can investigate the dynamical evaluation pro- 
cess using the deduced filter function: 

F{uL) = 8sin 4 ( — )sin 2 ( — )/cos 2 — (11) 
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Consider the frequency-dependent character of the 
evaluation process. From equation (8) it can be gen- 
erated that the elements of the resulting density matrix 
increase with a. Namely, decoherence is inhibited, in- 



stead of strengthen, by the optical frequency dispersion. 
This effect, though can be hardly observed due to the 
ordinary situation with a << wn, may more generally 
inspire us that under a DD control process an additional 
degree of freedom with special structure can stabilize the 
spin coherence. In the optical case here, since a signif- 
icant decoherence will emerge within a length scale of 
a 2 f(L) << 1, equation (8) can then be simplified as: 



Pl3, 14,23,24 = Pl3,14,23,24(0)e 



(12) 



which is not sensitive to the optical frequency dispersion. 

We now turn to the characterization of the degree 
of entanglement of the two-photons state. Wootters 
concurrence [24j is particularly convenient for the two- 
qubit case here. Other reliable measure of entanglement 
will yield the same conclusion. The concurrence can be 
calculated explicitly from the denity matrix p discribed 
in equation (12) 



C(p) — max(0, v A 



A 4 ) (13) 



where the quantities Xi are the eigenvalues in decreasing 
order of the matrix: 



(14) 



In the following we will analysis the evoluation of en- 
tanglement of a class of bipartite density matrices ini- 
tially prepared with the common X- form 25]: 

(pn pi4 N 

P22 P23 

P32 P33 

P41 p 4 4/ 

which occur in many contexts including pure Bell states 
as well as Werner mixed states. For the concurrence with 
such arbitrary initial state there is no compact analytical 
expression. However, it has been readily showed |26j that 
entanglement sudden death always occur under a clas- 
sic Gauss noise, as well as the evolution with DD control 
implemented. The DD method inhibit losing of entangle- 
ment by reshaping the exponential factor of the diagonal 
elements. We now give the numerical simulation to eval- 
uate the performance of the CPMG in our scheme. 

Without losing universality, we numerical analysis the 
concurrence evolution with a initial specific X-form en- 
tangled state: 

/1/2 

p (0) = 3 11 
\ 1/2/ 

Thus the initial concurrence is C(0) = -. In Fig. 2 
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(a) we give the numerical simulation of the concurrence 
evolution under a Gauss 1 /f noise (the noise spectral den- 
sity is S(uj) cx 1/w, origination of which are discussed in 
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FIG. 4: Concurrence evolution under a general Gauss noise 
environment with variable a ranges from 0.5 to 1.5 

Ref. [13] )• We observed that the concurrence of entangle- 
ment improves as the intervals between pulses decrease. 
As a contrast we also depicted the entanglement evo- 
lution process with SE pulse implemented in Fig. 2 (b). 
It can be inferred that the CPMG sequence drastically 
out-perform the single SE pulse in inhibiting the entan- 
glement sudden death. The performance of the SE pulse 
has been experimental tested as an effective method for 
coherence control! 281. 



We also analysis the entanglement preservation in a 
fixed fiber length as the total number of pulses changes. 
Thus the minimum number of pulses can be estimated to 
achieve a given level of entanglement. Fig. 3 shows that 
the entanglement rebirth while a certain number of pulses 
are implemented. However, it can be seen that the growth 
rate of concurrence declines while the number of pulses 
increases. This makes a preservation of entanglement to 
a extremely high level be difficult. 

The performance of CPMG sequence under a general 
Gauss l/f a with spectral density S(lu) oc l/uj a is also 
considered. Fig. 4 shows the concurrence evolution under 
the noise spectral density with variable a ranges from 0.5 
to 1.5. 

In conclusion, we demonstrated that dynamic decou- 
pling can effectively overcome the stochastic dispersion 
in the polarization-maintaining birefringent fibers under 
classic Gauss noises. Entanglement can be successful pre- 
served by embedding waveplates as CPMG sequence into 
the fibers. Our work will evidently enhance the scope of 
fiber-based quantum communication. We also hope this 
method be extended with more complicated DD sequence 
to general decoherent environment in ordinary sigle mode 
fibers. 
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